SHEP 01/01 



Lattice .B-PhysicsQ 



o 
o 



o 

-f— > 



Abstract 



C T Sachrajda 

Department of Physics and Astronomy 
University of Southampton 
O ! Southampton S015 7PW 

Sn' UK 
a 
m 

I review the status of lattice simulations relevant for phenomenological studies of 
£>. ■ 5-physics. Results for much-studied quantities such as fs, Bg and form- factors for 

semileptonic decays are presented as well as those for quantities which have begun to 
be studied only recently (such as S-lifetimes). The improvement in the precision of 
the determinations of the mass of the 6-quark, which has been made possible by new 
results for perturbative coefficients, is discussed. Finally I describe new ideas being 
developed with the aim of making computations of two-body non-leptonic decays 
possible. 



Jh : 1 Introduction 

> 

^ ■ At this conference we have seen a huge amount of experimental data being presented 
on 5-physics, data from which we would like to extract fundamental information about 
the parameters of the standard model and CP-violation, and (by overconstraining the 
comparison with the standard model) to search for signals of new physics. Of course much 
more data will become available in the coming years. A major difficulty in achieving the 
above goals is our inability to quantify non-perturbative strong interaction effects. Lattice 
QCD provides the opportunity for evaluating these effects ab initio, and there is indeed a 
wide-ranging program of numerical simulations of 5-physics. In this lecture I will briefly 
review some of the recent progress in lattice studies of I?-physics. 

It is not appropriate for me in this talk to discuss technical aspects of lattice calculations in 
any detail, I will therefore concentrate on presenting results and discussing the potential for 
future calculations. I start with some basic facts about lattice calculations in 5-physics. In 
any simulation the computing resources are limited, and one therefore has to compromise 
when minimising competing systematic errors, those due to the granularity of the lattice 



1 Invited lecture presented at the 7th International Conference on B-Physics at Hadron Machines, Beauty 
2000, Magaan, Sea of Galillee, Israel, 13-18 September, 2000. 
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(discretization effects) and finite-volume errors. In current simulations, this typically leads 
to a choice for the lattice spacing (a) in the range given by a -1 ~ 2-4GeV, i.e. a is 
larger than the Compton wavelength of the 6-quark. This means that we cannot study the 
propagation of a physical 6-quark directly and either have to i) use effective theories, such 
as the Heavy Quark Effective Theory (HQET) or Non-Relativistic QCD (NRQCD) or ii) 
calculate physical quantities with the heavy-quark mass tuq in the region of m c (the mass 
of the charm-quark), and perform the extrapolation to itlq = m^. 

The precision of lattice calculations is limited by systematic uncertainties. These can be 
overcome in principle with increased computing resources and much effort is being devoted 
to reducing these uncertainties (and in particular to eliminating quenching, i.e. to include 
vacuum polarisation effects). Most, but not all, of the results presented below have been 
obtained in the quenched approximation. It should be remembered however, that even in 
unquenched simulations the masses of the sea quarks are large |. Controlling the effects 
of sea-quarks is perhaps the major challenge for the lattice community and in recent years 
we have began to address this challenge. 

From lattice simulations we know how to calculate matrix elements of the form (H | O(0) | 0) 
and (Hi \ O(0) \ H 2 ), where H and Hi 2 are hadrons and O(0) is some local composite 
operator (several examples of such matrix elements will be discussed below). We are not 
yet in a position to compute matrix elements in which there are two or more hadrons in the 
initial or final state, and are therefore unable to compute B — > H\Hi decays amplitudes, 
but progress is being made in this area (see sect. [5] below). 

The plan of this lecture is as follows. In the next section I briefly review the results for 
some standard physical quantities: leptonic decay constants fs and js a \ the ^-parameters 
of B-B mixing and form-factors for semileptonic decays. Sect. contains a discussion 
of computations of inclusive decays and lifetimes. Results have been recently obtained 
for these quantities but they have not been studied as extensively as those in sect. ^|. 
New calculations of perturbative coefficients have enabled the mass of the 6-quark to be 
determined from lattice simulations with considerably better precision and this is reviewed 
in sect. f|. I then describe recent attempts to develop techniques to enable two-body 
non-leptonic decays to be studied on the lattice (sect. |^). Finally I summarise the main 
conclusions in sect. |[ 

2 Decay Constants, B— Parameters and Form-Factors 

Leptonic decay constants J'b and fs s , the .B-parameters of B-B mixing and the form- 
factors for semileptonic -B-decays have been computed in lattice simulations for over ten 

2 For example we have yet to simulate with sufficiently light quarks so that the p — > hit decay is 
kinematically possible. 
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Figure 1: Diagrammatic representation of a leptonic decay of a S-meson. 



years now. In this section I briefly discuss each of these quantities in turn. 



2.1 f B 

The strong interaction effects in leptonic decays of 5-mesons (see fig.[l]) are contained in 
the matrix element 

(0\A^0)\B(p)) = tf BPfl . (1) 

Lorentz and Parity Invariance imply that all the non-perturbative QCD effects are para- 
metrized in terms of a single number, f B , defined in eq.(|TJ). 

Claude Bernard, the reviewer at the Lattice 2000 symposium |l] summarised the current 
status of the results as Q: 

Full Theory N f = 



f B 200 ± 30 MeV 175 ± 20 MeV 

fsJfB 1.16 ±0.04 1.15 ±0.04 



At last year's conference the reviewer S.Hashimoto stated that all available data is consis- 
tent with the following estimates: 0] 

Nf = 2 N f = 



f B 210 ± 30 MeV 170 ± 20 MeV 

f Bs 245 ± 30 MeV 195 ± 20 MeV 
f B Jf B 1.16 ±0.04 1.15 ±0.04 



These results have been largely stable for some time. The last time that I performed a 
compilation was in 1997 together with Jonathan Flynn || and reported (on what were 
mainly quenched calculations): 

f B = 170 ±35 MeV, f Bs = 195 ± 35 MeV, ^ = 1.14 ±0.08. (2) 

JB 



3 Ref. Q contains a detailed critical analysis of all the lattice results for the quantities considered in 
this section. 
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Although the results have been stable, the errors have been decreasing very slowly (if at 
all). The errors will not decrease significantly until we begin to get a serious control of 
quenching errors. In unquenched calculations, the value of the lattice spacing typically 
varies by 10% or so depending on the physical quantity which is used to set the scale. It is 
therefore not really possible to determine dimensionful quantities such as /# with a better 
precision than about 10%. 

In the last two years we have began to study quenching effects, although it must be 
remembered that at present the masses of the sea-quarks are still relatively heavy (typically 
a little smaller that the strange quark mass). These effects can only be studied meaningfully 
if all other variables are kept constant. Using results from the MILC ||], CPPACS |5j and 
NRQCD || collaborations, C. Bernard estimates that there is about a 20 MeV increase in 
the value of fs in going from the quenched approximation to the Nf = 2 case (using f n to 
set the scale) 0. The value of /b s //b appears to be the same in the two cases. 

Although the errors have not been reduced substantially, there have been many systematic 
checks on the stability of the results. For example discretization errors have been studied 
by using improved actions and/or extrapolating the results to the continuum limit. It 
had been suggested that may be significantly lower than the results presented above 
because of discretization errors (this suggestion was based on some extrapolations to the 
chiral limit), but this has not survived more careful analyses |L|. 

When lattice computations of the leptonic decay constants of heavy mesons were beginning 
there were no experimental data. In 1997 we compiled the (quenched) lattice predictions 
for the decays constants of charmed mesons as ||: 



The Particle Data Group's review this year quotes f Ds = 280 ± 19 ± 28 ± 34 MeV as 
the best experimental result (and recent results, still with large errors, are compatible with 
this). It will be interesting to observe future lattice and experimental results as the errors 
in both decrease. 

Finally let me mention that lattice computations have shown that there are large correc- 
tions to the HQET scaling law for the decay constants of heavy-light pseudoscalar mesons P 



f D = 200 ± 30 MeV and f Ds = 220 ± 30 MeV. 



(3) 



fp 



1 



(4) 



2.2 B B : 



The relevant matrix element for B -B mixing is: 



M(ii) = ( B° I 67^(1 - 7s)g W - 7s)g \B°), 
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Figure 2: Diagramatic representation of B -B mixing, 
where q represents d or s. This is illustrated in fig. |2|. 

Following the conventions introduced in kaon physics, B-parameters are defined by 

M{fi) = \fW B B B ^). (5) 

B B {jj) is scheme and scale-dependent, so it is convenient to define a renormalization scheme 
and scale- independent (up to NLO) quantity 



Bt = a s (fi) 2 ^ 



B B (ji). (6) 



where J n is a known constant. Many groups have studied mixing in lattice simulations 
and C. Bernard, the reviewer at this year's lattice conference, summarised the results as: 

B Bd = 1.30 ±0.12 ±0.13, B B jB Bd = 1.00 ±0.04, (7) 



f B \fi~ B ~ d = 230 ±40 MeV, g= B ° = 1.16 ± 0.05 . (8) 

The second error for B Bd is an estimate of the error due to quenching. £ is a particularly 
important quantity for studies of the unitarity triangle. It should be noted that the error 
on £ is not a small one from the lattice perspective, indeed given its central importance in 
phenomenological studies most groups tend to be cautious in determining this error. The 
reason that it should not be considered small is that the quantity which is being computed 
is £ — 1 = 0.16 ±0.05, i.e. it has a quoted error of over 30%. In my judgement therefore the 
uncertainty in the value of £ should not be further inflated in phenomenological studies. 
The results for £ are stable, and there is no evidence that they change as sea-quarks are 
introduced. 

For comparison with the results quoted above I present also those quoted in our 1997 
review |§: 

B Bd = 1.4 ±0.1, B Bd /B B , = lM±0M, (9) 



— Sb \ b b 

f B \jB Bd = 201(42) MeV, g= ' v ' = 1.14 ± 0.08. (10) 

f B y B B 
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These results have been obtained with propagating heavy quarks (i.e. by extrapolating 
results obtained with mq ~ m c ). In the static approximation the perturbative corrections 
are very large and hence these have not given us much useful information to date. 

2.3 Exclusive Semi-Leptonic .B-Decays 

Exclusive semileptonic decays of 5-mesons are being used to determine the V c b and V u \, 
CKM-matrix elements. Diagramatically the amplitudes can be represented by the diagram 
in fig. |||. Lorentz and parity invariance imply that the amplitudes can be expressed in terms 




Figure 3: Diagramatic Representation of Semileptonic B-Decays. 



of invariant form factors. For example for the decays into a pseudoscalar P (B — > 7r or 
B — > D decays), in the helicity basis, 

(P( PP ) |V M (0) I B(p B ) ) = /V) Ml ' q2 Ml q, + f + (q 2 ) {(pb+Pp), - Ml ~ 2 M2p %) ■ 

(11) 

As a result of parity invariance only the vector component (from the V — A weak current) 
contributes in this case. When the final-state hadron is a vector meson (p or D*) the 
amplitude is written in terms of four invariant form factors. 

I will only make a few brief comments about semileptonic decays: 

• Lattice simulations of B — > n, p decays require the momentum of the final-state light 
meson to be small in order to avoid discretization errors. This means that from 
lattice simulations of B — > 7r,p semileptonic decays we only obtain results at large 
values of q 2 . This is illustrated by the first diagram in fig. (|, which contains recent 
data from the APE [| and UKQCD collaborations 0. 

• Lattice computations of semileptonic form-factors have been performed for many 
years now. At the lattice 2000 symposium updated results were presented from the 



The curves are APE collaborations fits to their points 
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Figure 4: The left-hand diagram shows results for /+ (upper points) and f° (lower points) 
from the APE || and UKQCD || collaborations. The curves are fits to the points from 
APE, showing an extrapolation to lower values of q 2 . The right-hand diagram shows 
the UKQCD results for the distribution for b — > p decays in the range 14GeV 2 < q 2 < 
20.3 GeV 2 , which corresponds to the high q 2 bin of CLEO. The region marked phase space 
only is inaccessible to the lattice calculations and the curve corresponds to taking the 
form-factors from the last lattice point and extrapolating using the phase-space only. 



APE, UKQCD, JLQCD and Fermilab collaborations and "the results were fairly 
consistent in the region where all groups had direct calculations (19 GeV 2 < q 2 < 
23 GeV 2 )" 0. As an example, the first dia gram in fig. [| shows the results for the 
form factors / + and f° for B — > 7r decays. 

Much effort is being devoted to extrapolating the lattice results to smaller values of 
q 2 , using as many theoretical constraints as possible (e.g. HQET scaling relations, 
unitarity and analyticity, kinematical constraints, soft pion relations etc.) flj~0"l . The 
results of such extrapolations for f° and / + are also shown in fig.^. 

For a direct application of lattice computations one should compare the lattice results 
to experimental distributions at large values of q 2 . An example of such a comparison 
is shown in the second diagram in fig.[|, where the lattice results from the UKQCD 
collaboration for B — > p decays are compared with the high-g 2 bin from CLEO jTTJ . 
Integrating the results in this bin the comparison takes the form: 

AT(14 < g < 20.3) = UKQCD preliminary (12) 

ps 1 GeV 

= 7.1(2.4) x 1(T 5 CLEO (13) 
from which one obtains V u b = 2.9(0.5) x 10~ 3 . 
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3 Inclusive Decays and Lifetimes of Beauty Hadrons 



In this section I will briefly mention some physical quantities for which the first lattice 
calculations were performed in the last few years. It will certainly be possible to improve 
on the precision of these first-generation computations. 



3.1 Lifetimes of Beauty Hadrons: 

The fact that the 6-quark is heavy makes it possible to derive an operator product expansion 
for inclusive beauty decays [|l2j, which results in an expansion in inverse powers of m b . 



Specifically for widths of beauty hadrons: 

r ™ = 192.3 (H + ) wJ I ( } 

where c 3 and c 5 are calculable in perturbation theory and 

-Ai = — !— (H b \hD 2 h\H b ) and 3A 2 = — ?— ( F 6 1 f ha^h \ H b ) (15) 

defined in the rest frame of H b (Aj and A 2 are the matrix elements of the kinetic energy 
and chromomagnetic operators respectively). This leads to 

T ' B ' ] l + 0(>) (16) 



t(B°) \m 

r(A b ) XMb) ~ *i(B°) A 2 (A b )-A 2 (£ 

r(B°) 2m 2 b m l 



+ 3c G zv 0/ 2 zv ' +••• (17) 



(0.98 ± 0.01) + O ^ , (18) 



to be compared to the experimental values f\ 

r(B~) r(A 6 ) , , 

v ; = 1,06 ±0.04 and = 0.79 ± 0.05 . (19) 

T[B") T[B") 

In view of the discrepancy between the theoretical prediction in (|18D and the experimental 
results for the ratio r(A b ) /r(B°) it is important to evaluate the 0(l/m\) corrections 0. 
These may be significant because it is only at this order that the light-quark in 5-meson 
is involved in the weak decay, and hence their contribution is enhanced by phase-space 

5 See also ref@ in which the result t{B-)/t{B°) = 1.07(2) is quoted. 

6 The significance of the discrepancy is underlined when one notes that it is really 1 — t(A{,)/t(B°) 
which we calculate. 
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factors [H|. To evaluate these 0(1/ mf) corrections we need to compute the matrix elements 



of the following dimension-6 operators ||14| : 

Oi = h,(i-i 5 )q qi"(i-i 5 )b, o 2 = &(i- 7 5 )g ?(i + r% 

T l = b 1 ^-l 5 )T a qq^{l-l 5 )T a b and T 2 = 6(l- 7 5 )T a g q(l + 7 5 )T a 6, (20) 
where T a represents the colour matrix. 

For mesons, the evaluation of these matrix elements is very similar to that of B B and is 
totally straightforward. For the A b baryon the calculation is a little less straightforward, 



but can be performed using standard techniques. Together with M. di Pierro we found ||T5 

' lL ' -- 1.03 ±0.02 ±0.03 , (21) 



r(B d 



in good agreement with the experimental result. The first error in fl2T|) is the lattice one 
and the second is our estimate of the uncertainty due to the fact that the (perturbative) 
coefficient function is only known to one-loop. There is a remarkable factorization of the 
lattice operators, i.e. the lattice 5-parameters are close to 1. For the A b baryon we have 
only performed an exploratory calculation, in which we did not have sufficient control of 
the behaviour with the mass of the light-quark [|l6j. The spectator effects were found to 



be very significant but, at least for lights quarks corresponding to m n ~1 GeV, the result 
of r(A b ) /t(Bq) — 0.91 would imply that spectator effects would not be sufficiently large 
to account for the full discrepancy. It is clear however, that these calculations need to be 
improved (which will be relatively simple to do) and that this will make a considerable 
impact on our understanding of this important question. 



3.2 £? s -Meson Lifetime Differences: 

We now consider the difference of the widths of the B s mesons. Using the operator product 



expansion we have |L7 



Ar, Glm 2 



V cb V cs \ T Bs x 



T Bs 127im Bs 

{G(z)(B s \Q L (m b )\B s )-G s (z)(B s \Q s (m b )\B s )+5 1/m VT^} (22) 

where z = m c /m b , the operators are Ql = bq^s bj^s and Qs = b(l — 7 s)s 6(1 — 7 s)s, 
the 5i/ m term represents 0{l/m b ) corrections and G(z) and Gs(z) have been computed in 
perturbation theory up to NLO [Q. The evaluation of these hadronic matrix elements has 



recently began to be performed in lattice simulations by Hashimoto et al. [19j and APE 
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collaboration |2D N. The APE collaboration quote 



= (4.7 ±1.5 ±1.6)% APE 0, (23) 

r - 



where the second error represents an assumed 30% uncertainty on the 1/rrib corrections. 
The calculations were obtained with propagating heavy quarks. The result from Hashimoto 
et al. is 



Ar 



B. 



B s 



;i0.7±2.6±1.4±1.7)% Hashimoto et al. @, (24) 



where the first error is due to the uncertainty in our knowledge of fg a , the second is the 
lattice error in the matrix elements and the third is the estimate of the error due to the 
0(l/rrib) corrections. These results were obtained by simulating NRQCD. 



The central values in eqs. ( p3|) and (24) are somewhat different (although they are compat- 
ible within errors), however it should be noted that this difference is not due to different 
lattice results for the matrix elements. Indeed the values of the matrix elements deter- 
mined by the two groups are in good agreement. The difference is due to the different 
inputs: Hashimoto et al. use a lattice result for fg s = 245 ± 30 MeV and the experimental 
value of the inclusive semileptonic branching ratio; APE use the lattice value of £ and the 
experimental value of r Bs Arne d mB a /mB d - 



4 The mass of the b- quark 

The mass of the 6-quark, nib, is one of the parameters of the standard model. In this 
section I will review recent progress in evaluating rrib from simulations of effective theories, 
such as the HQET or NRQCD Q. The discussion will be based on the evaluation of the 
two-point correlation function 

C(t) = (0\A o (x,t)A Q (0,0)\0) (25) 

in the HQET. is the axial current, A^ = h^^q, and h (q) is the heavy-quark (light- 
quark) field. The evaluation of C(t) is relatively straightforward and has been performed for 
many years now. The progress in recent years has been in the theoretical understanding of 
how one can determine mj, from C (t) , and in the perturbative calculations which have been 
performed up to three-loop order enabling the results to be obtained with good precision. 

7 There are also some even more recent results in the static limit, both in the quenched approximation 
and with two flavours of sea quarks ^T|. The results quoted in these papers are 5. 1±1.9±1. 7% (quenched) 
and 4.3 ± 2.0 ± 1.9% (n f = 2), to be compared to those in eqs.Q and (§4|). 

8 For a detailed review of lattice determinations of quark masses in general, and mj in particular, see 
ref. @. 
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At large times t 



C(t) ~ Z 2 exp(-£t) (26) 

and from the prefactor Z we obtain the value of the decay constant /# in the static 
approximation f2~3fl . It is for this reason that calculations of C(t) have been performed for 



some years now. It was realised later that one can obtain m b from the measured value of 
£ (up to 0(A^ CD /m b ) corrections) |2Tf . 

The relation between £ and m& is a delicate one: 

Mb = m P b ° lc + £-5m, (27) 

where M# ~ 5.28 GeV is the mass of the .B-meson, m? olc is the pole-mass of the 6-quark 



and £ is determined from lattice computations using eq. (|26|) . 5m is the residual mass 



generated in perturbation theory in the lattice formulation of the HQET, 

8m = ~Y J XnOt n s {m b ) . (28) 

n 

In other words even if the action in the HQET is written without a mass term, higher order 
terms in perturbation theory generate such a term, 5m. Each term in the series for 5m 
diverges linearly as a — > 0, and these terms (partially) cancel the divergence in £. In order 
for the cancellation to be sufficiently precise the lattice spacing a cannot be too small. 

The perturbation series for 5m diverges. Indeed it is not Borel Summable, leading to 
a renormalon ambiguity, which is an intrinsic ambiguity of O(Aqqtj). The pole mass 



of a quark is also not a physical quantity; it also contains a renormalon ambiguity [25 



26]. The renormalons in the pole mass and 5m cancel (for a detailed discussion of the 
cancellation of such renormalons see ref. j?/]). Let m b be the 6-quark mass defined in the 
MS renormalization scheme defined at the scale mj itself (I use m b as an example of a 
"physical" quark mass). The pole mass and m;, can be related in perturbation theory: 

m b = [ 1 + J2 Dna ^ mb ^ m b° le ( 29 ) 

n 

= [1 + Y,D n a n s {m b )\ {M B -£ + ^X n <(m 6 )} . (30) 



4.1 Perturbative Matching 

All the perturbative coefficients necessary to compute to N 3 LO have recently been 
calculated. 

• The coefficients D 2 and D 3 in the relation ( p9|) between rn b and the pole mass have 
been calculated in refs. ]28| and respectively. 
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The relation between the lattice (bare) coupling constant «o and the MS coupling 
takes the form 

a s (/i) = 1 + ^2d n (a/j)a>Q a . (31) 



and the coefficient d 2 can be found in refs |31| . 

Finally we need to consider the perturbative expansion of the residual mass: 

aSm = Citt + C 2 al + C 3 al H (32) 

— C\ is well known, C\ = 2.1173. 

— C 2 has been determined from the calculation of large Wilson loops W(R, T), for 
which 

W(R, T) ~ exp(-25m (R + T)), (33) 

where R and T are the spatial and temporal lengths of the Wilson loop. 

C 2 = 11.152+n / (-0.282+0.035c sw -0.39l4w) H, where c sw is the coefficient 



of the Sheikholeslami-Wohlert term in the lattice action 133 



The Parma-Milan group has pioneered the numerical computation of perturba- 
tive coefficients at high orders using the stochastic formulation of theory. For 
the coefficients in eq. (|32|) they find K 



d = 2.09(4), C 2 = 10.7(7), C 3 = 86.2(5) . (34) 

From a numerical Monte-Carlo study on very small and very fine-grained lattices 
Lepage at al. |35|] find: 



C 3 = 80(2) + 5.6(1.8) . (35) 
For NRQCD only the first coefficient (Ci) is known. 



4.2 Numerical results for 



Before presenting the numerical results I repeat that the evaluation of £ is relatively 
straightforward. The recent progress has been in the evaluation of the high order per- 



turbative coefficients described in subsection |4.1| . The results obtained using the method 
described in this section have been: 



4.15 ± 0.05 ± 0.20 GeV |6 



(36) 



obtained in quenched QCD at NLO (i.e. with only the one-loop coefficient C\). 
The second error is the estimate of the uncertainty due to higher order perturbation 
theory. 
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m b = 4.30 ±0.05 ±0.10 GeV ||, (37) 
obtained in quenched QCD at N 2 LO (with coefficients up to Cq). 

m b = 4.30 ± 0.05 ± 0.05 GeV |J7|, (38) 
obtained in quenched QCD at N 3 LO (using all the known coefficients of section |4~T 



m b = 4.34 ± 0.03 ± 0.06 GeV |8 



obtained from the static limit of quenched NRQCD at N 3 LO. 
There is also an unquenched result: 

m b = 4.26 ± 0.06 ± 0.07 GeV g9[, 



obtained in unquenched QCD, with 2 flavours of sea-quarks at N 2 LO. 



(39) 



(40) 



We take the results in eqs. (|38|) (or ([39])) and (fiTf ) as the current best estimates for m b . 



Results obtained using NRQCD are in good agreement with those in the static theory, but 
we need to understand the errors due to higher orders of perturbation theory in that case. 



5 Nonleptonic Decays 



At this meeting many interesting results have been presented for exclusive two-body decays 
of S-mesons. Indeed these decays are one of the principal sources of information about the 
unitarity triangle and CP- violation. Apart from the golden process B — > J/^K S , the de- 
termination of the fundamental parameters from two-body 5-decays is severely restricted 
by our inability to quantify the non-perturbative strong interaction effects. At present we 
are some way from understanding how to compute these effects in lattice simulations, how- 
ever considerable effort is being devoted to developing the necessary theoretical framework. 
Progress has been made in understanding how to compute K — > -nit decay amplitudes 




and I will briefly discuss recent ideas 
need to be understood: 



41, 42|. There are two important features which 
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1. Lattice calculations are performed in Euclidean space and hence yield real quantities. 
One can therefore ask how one can get the full decay amplitude including the phase 
due to the final state interactions? It is true that from each correlation function 
one obtains a real quantity, however different correlators give different quantities and 
the full amplitude can (at least in principle) be reconstructed. For example, one 
can obtain the modulus of the amplitude from one correlator and the real part from 
another which is sufficient to determine the amplitude |42| . 



2. Although momentum is conserved in lattice calculations of correlation functions, 
energy is not (correlation functions are not integrated over time). 

p = 

K 




In this diagram the kaon decays at the origin producing two pions, which are annihi- 
lated at t\ and t 2 respectively with momenta q and —q. However the pions interact, 
and can rescatter (as represented by the grey oval). Whilst the total momentum 
remains zero, energy is not conserved and hence there is a contribution in which the 
two pions are each at rest (this is the lowest energy state). Correlation functions 
decay exponentially with time as exp(-Et), where E is the energy, and hence at 
large times the dominant contribution comes from the unphysical matrix element 
corresponding to a kaon at rest decaying into two pions, each of which is at rest. 
This is a major difficulty for lattice calculations. 

In a finite volume the energy levels of the two-pion states are discrete and Luscher 
and Lellouch have proposed to exploit this fact to determine the decay amplitude 
directly j|T|. One needs to have a state for which the energy of the two-pion state is 



equal to tuk, the mass of the kaon, and to be able to isolate the term corresponding 
to this energy. In practice, for the near future, this is likely to be the first excited 
state. This would require a lattice of about 6fm, which is large but not hugely so (]. 



The ideas discussed above need to be explored numerically and developed further. In par- 
ticular we need to understand how to include the contributions from inelastic channels 
(i.e. contributions from intermediate states other than two-pion ones). The current nu- 
merical studies of K —>■ mi decays are performed on lattices which are considerably smaller 
than 6fm, and so correspond to unphysical kinematics. However, one can then use chi- 
ral perturbation theory to obtain the physical amplitudes. For 5-physics of course chiral 
perturbation theory is inapplicable and the neglect of inelastic channels is not possible, so 
more theoretical developments are necessary. 

9 The finite volume corrections have also been discussed in considerable detail in refs. Ell, fr^l. 
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6 Conclusions 



I conclude with a brief summary of the main points of this talk: 

• Lattice simulations are a central tool in the determination of information about 
fundamental physics from experimental data. They are being applied to a wide 
range of quantities in 5-physics. 

• The numerical results presented in sections 0-|4] are an important part of the conclu- 
sions. 

• For standard quantities, such as fs and B B , the emphasis is on the control of sys- 
tematic errors, particularly those due to quenching. For other quantities (e.g. the 
applications to inclusive decays discussed in section ||) the calculations are just be- 
ginning. 

• For B — > light-meson semileptonic (and radiative) decays, the matrix elements are 
computed at large values of q 2 and theoretically motivated extrapolations are per- 
formed to lower values of q 2 . Where possible, to make optimum use of lattice results, 
a comparison with the experimental distributions at large values of q 2 should be 
made. 

• For two-body non-leptonic decays theoretical developments are needed before lattice 
simulations can provide information. Progress has been made for K —>■ tttt decays 
and we can look forward to extensions to 5-physics. 
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